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A THEORETICAL ANALYSIS OF THE EFFECTS OF FUEL MOTION ON AIRPLANE DYNAMICS ' 


By ALBERT A. Segv 


SUMMARY 


The general equations of motion for an airplane with a num- 
ber of spherical fuel tanks are derived. The motion of the fuel is 
approximated by the motion of solid pendulums. The same type 
of derivation and equations are shown to apply to any type of 
fuel tank where the motion of the fuel may be represented in 
terms of undamped harmonic oscillators. 

-Motions are calculated for a present-day high-speed airplane 
and a free-flying airplane model with two spherical tanks in the 
symmetry plane. These calculations show that the normal air- 
plane motion may be considerably modified and that residual 
oscillations may result. The ratio of the natural fuel fre- 
quency to the natural airplane frequency is shown to be the most 
important parameter for determining the effect of the fuel motion 
on the airplane motion. The stabilizing effect of turbulence in 
the fuel is discussed, and it is suggested that the stabilizing effect 
of artificially induced turbulence be investigated experimentally. 


INTRODUCTION 


Small-amplitude lightly damped lateral oscillations are a 
troublesome characteristic of certain high-speed airplanes. 
Several possible explanations for these oscillations, which are 
adequate in specific cases, have been offered. For example, 
reference 1 shows that nonlinear aerodynamic derivatives 
could cause such oscillations, and it has been shown that 
atmosphefic turbulence is another possible cause. It has 
also been suggested that a possible cause of such oscillations 
is the motion of fuel in the tanks. In some recently designed 
airplanes the mass of the fuel relative to the airplane mass is 
much larger than has been common in the past; therefore, 
the effects of fuel motion ean be expected to be relatively 
more important. In fact, in several cases baffling the fuel 
tanks was found to have considerable effect on the general 
handling qualities of the airplane and sometimes actually 
eliminated the troublesome lightly damped lateral oscilla- 
tions which had been present. 

An experimental investigation of the effects of fuel motion 
on the lateral motion of a free-flying airplane model is de- 
scribed in reference 2. The results indicated that the effects 
of fuel motion were noticeable and caused the lateral motion 
of the model to be very erratic. 





The present analysis treats each fuel tank as a pendulum 
oscillating in two degrees of freedom and applies Lagrange's 
equations of motion to obtain the interaction between these 
pendulums and the airplane. Thus, for small motions the 
fuels are treated as simple harmonic oscillators. The results 
are applied to obtain the general equations of motion of this 
system and, in particular, the lateral motion of an airplane 
with internal fuel tanks in the plane of symmetry of the air- 
plane. Since the general solution of the equations is ex- 
tremely complicated, an attempt is made to evaluste the 
results by carrying out numerieal calculations for specific 
cases. This approach is shown to be adequate in yielding 
the most general effects of fuel motion. 

The discussion of the numerical application of the equa- 
tions of motion to specifie cases is given in detail after the 
derivation of the equations of motion. This discussion of 
results is understandable quite independently of the deriva- 
tion of the equations of motion. 


SYMBOLS 

X,Y,Z airplane stability axes with origin deter- 
mined by equations (13); also compo- 
nents of applied forces along these axes 

L, M, N components of applied moments about 
X-, Y-, and Z-axes, respectively 

i,j,k unit vectors along .X-, J-, and Z-axes, 
respectively 

T, Y, 2 components of translational displacement 
of airplane 

F vector translational velocity of airplane 
(ix + jy + kz iU v) 

U magnitude of steady-state velocity 

U, &, W components of disturbance translational 
velocity of airplane 

voo vector disturbance velocity of airplane 
(iu-+jo+ kw) 

R vector position of a point in airplane 
GR. T jR,T KR) 

Ry vector position of center of gravity of 
fuel in a particular fuel tank 

y total vector velocity of a point in airplane 


‘Supersedes NACA TN 2280, “A Theoretical Analysis of the Effects of Fuel Motion on Airplane Dynamics” by Albert A. Schy, 1951. 
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aideslip angle (tan v) 


infinitesimal rotations of airplane about 
X-, Y-, and Z-axes, respectively 

vector rotational velocity of airplane 
(ig + ó-- ky) 

components of angular displacement from 
vertical of line joining fuel center of 
gravity to tank center, taken on mutu- 
ally perpendicular planes; ¢ positive in 
direction of positive roll and 7 positive 
in direction of positive pitch 

distance from tank center to fuel center 
of gravity 

vertical displacement of fuel center of 
gravity from equilibrium position 

number of fuel tanks 

mass 

total mass of airplane and fuel 

total moments and produet of inertia of 
airplane about X-, Y-, and Z-axes 

rigid-body moments and products of in- 
ertia about axes through center of gravity 

fuel moments of inertia about f- and 
y-axes through tank center 


fuel moments of inertia about ¢- and 
y-axes through fuel center of gravity 
nondimensional radius of gyration in roll 


(és 


nondimensional radius of gyration in yaw 


( ae 
mb? 


nondimensional product-of-inertia param- 


y 
eter CE) 


kinetic energy 

potential energy 

period 

time for exponentially damped or increas- 





ing oscillation to halve or double am- 


plitude, respectively 
time 


; 2 ; "Ut 
nondimensional time parameter (=) 


acceleration due to gravity 


flight-path angle with respect to horizontal 
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air density 

wing area 

wing span 

lateral nondimensional mass coefficient 


( m, mg 
EP Moy Sb 
differentiation operator ds) 


trim lift coefficient map eese l 
3? U*S 
rolling-moment coefficient 


Rolling moment 
} pU'Sb ) 


yawing-moment coefficient 


Yawing moment 
5 oU*Sb ) 


] Lateral force 
lateral-force coefficient ir 


l Us 


particular fuel tank, or summation index 
over fuel tanks (f=1, 2, . . . k) 

airplane without fuel 

initial conditions at £=0 

















A THEORETICAL ANALYSIS OF THE EFFECTS OF FUEL MOTION ON AIRPLANE DYNAMICS 


DERIVATION OF EQUATIONS OF MOTION 
ASSUMPTIONS FOR DERIVATION OF GENERAL EQUATIONS OF MOTION 


As a first approximation, only the effect of the motion of 
the fuel as a whole is considered; that is. only the funda- 
mental mode of the wave motion is considered. and this 
mode is approximated by rigid-body motion. The main 
effect of the internal wave,motion is to introduce damping 
into the fuel oscillation. This damping is caused by the 
conversion of kinetie energy into heat through the turbu- 
lence caused by the splashing of the fuel. A strictly analytic 
consideration of such damping effects is extremely difficult; 
on the other hand. the damping caused by the viscous tan- 
gential forces between the fuel and the tank is completely 
negligible (see reference 3). The analysis of the problem is 
therefore confined to the motion with no fuel damping and 
the effect of the damping is considered in the discussion of 
the results. 

In a spherical tank the fuel can oscillate approximately as 
a rigid body if no splashing is assumed for small oscillations. 
The motion may be pictured as the “rocking” of a spherical 
segment of constant shape. The restraining force of the 
tank, which always acts in a direction normal to the motion, 
is exactly analogous to the tension in a pendulum. Thus, 
the small motions of the fuel in a spherical tank may be rep- 
resented quite well by the well-known simple properties of 
small pendulum motions. This approach is used in the 
mathematical analvsis of the problem. 

The effect of aspherical tank shape can be &pproximated 
by replacing the tank by an equivalent harmonie oscillater 
with an arbitrary amount of turbulence damping added even 
for small motions. For example, the representation of rec- 
tangular tanks as harmonic oscillators is discussed in refer- 
ence 4, Thus in this ease also the most general effects of 
the fuel motion on the airplane motion should be qualita- 
tively obtainable by this type of analysis. 

The effects of large-amplitude fuel motions will be dis- 
eussed qualitatively after the discussion of the results of the 
mathematical analysis. As usual in stability analysis ell 
motions are assumed small and second-order terms are 
ignored, 


DERIVATION OF GENERAL EQUATIONS OF MOTION 


With the preceding assumptions the physical problem can 
be considered as the interaction between two or more rigid 
bodies, namely the airplane and the several fuel pendulums, 
with each fuel pendulum considered as suspended from the 
tank center. The only potential energy considered in the 
system is that of the pendulums. If the inertial character- 
istics of the airplane and the fuel are known, the kinetic 
energy of the system can be obtained from the translational 
and rotational velocities of the airplane and the fuels. With 
this information the interactions in the system can beobtained 
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by using Lagrange's equations of motion in the form (see 
reference 5) 


d (DEA OE, OE, 

5 (SA) x-- 2B 9; o 
dt NOg, Òg: 9g: 

where q; is one of the n generalized coordinates of the system 
corresponding to the n degrees of freedom, @, is the corre- 
sponding velocity, and Q, is the corresponding generalized 
force. The g, will be lengths and angles and the correspond- 
ing Q; will be forees and moments, respectively. 

The airplane itself introduces the customary six degrees 
of freedom, whieh are the three displacements of the airplane 
system along axes fixed in the airplane (z, y, z) and the cor- 
responding angles of rotation of the airplane about these 
axes (ó, 0, V). For small displacements, the pendulum 
motion can be described by two angles ¢ and y since the 
vertical motion can be neglected (see fig. 1). The angle ¢ is 
measured from a vertieal line through the tank center to the 
projection of the line joining the tank center to the fuel 
center of gravity on the vertical plane parallel to the Y-axis 
and ; is the corresponding angle in the vertical plane parallel 
to the X-axis. For small angles, t and 5 may be represented 
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FIGURE 1.—Decomposition of horizontal fuel motion in terms of angles ¢ and s» in vertical 
Planes. 
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as in figure 1. In effect this figure makes use of the fact that 
small angles may be added vectorially. When the two addi- 
tional coordinates ¢ and 7 are used to describe the pendulum 
motion, the whole system has two additional degrees of 
freedom for each fuel tank. 

Expressions must be obtained for E, and E, in terms of 
the coordinates of the system and their time derivatives in 
order to use equation (1). The only potential energy is that 
of the fuel pendulums, which can be written as follows for k 
fuel tanks: 


k 
E=}, mjh,g (2) 


For the height of the center of gravity in each fuel tank, as 
can be seen from figure 1, 


h=hyt+h,~I—I cos focos vet[1-(1 -i e) »)] 


or 


h elus) (8) 


Note that the vertical displacement A is of second order in 
the small quantities n and ¢. This fact justifies the previous 
statement that the vertical displacement could be neglected 
in describing the pendulum motion only by the two coor- 
dinates 7 and ¢. As might be expected, equations (2) and 
(3) indicate that each fuel pendulum is being considered as 
an undamped oscillator with two degrees of freedom in a 
horizontal plane. 

The kinetic energy of the total system can be written as 
the sum of the kinetic energy of the airplane and the kinetic 
energies of the fuels. Also the kinetic energy of each rigid 
body can be expressed as the sum of the translational energy 
of the mass moving with the velocity of its center of gravity 
and the rotational kinetic energy of the mass about its center 
of gravity. Thus, when the inertial characteristics of the air- 
plane and the fuels are known, the kinetic energy can be 
obtained as a function of the generalized coordinates and 
velocities if the translational velocity of each center of gravity 
and the angular velocities of the airplane and fuels about 
their respective centers of gravity can be expressed in terms 
of these generalized coordinates and velocities. 

In order to obtain the required expressions for these veloci- 
ties, a system of axes fixed in the airplane with the X-axis 
along the steady-state velocity at t=0 is used, asis customary 
in stability analysis. For the present the origin of the coor- 
dinates will not be specified. However, these stability axes 
are not inertial axes and Newton’s second law applies only in 
ən inertial system of axes. The inertial axes may be taken 
as axes fixed in the earth. Then in the equations of motion 
the velocities and accelerations must be measured with Te- 
spect to the earth, and their expressions in terms of com- 
ponents in the moving airplane axes may be obtained as 
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shown in reference 6. These expressions will give the kinetie- 
reaction forces, which for the case of a rotating system are 
often referred to as “gyroscopic” forces. For the velocity, 
referred to the inertial system, of any point defined by the 
vector R in the airplane axes (in particular, for the centers 
of gravity previously discussed), 


V=F+oXR+R=iU +0+0XR+R (4) 
where all vectors are given in terms of the airplane axes and 
Feit+jytki=iU +o 


is the velocity of the origin of the airplane system with respect 
to the earth, while v and w are the translational and rota- 
tional disturbance velocities of the airplane axes. 

Equation (4) may now be used to express the inertial 
velocities of the airplane center of gravity and the fuel 
centers of gravity in terms of the generalized coordinates by 
inserting for R the values Ra and R, where Rais the vector 
position of the airplane center of gravity and R/ indicates 
the vector position of any particular fuel center of gravity. 
The vector R, is constant; therefore, R,—0. To obtain È’, 
note that to first order 


R/ =R,— il, sin (ysd-6— 1) — jitro cos yo— y sin yot 
kl, cos (yot @— 1) 


R/~R,—il,sin Yot(@— ns) cos qd —FIASs— COS Yo— 
y sin y)-+kl,[cos yo— (8— ny) sin Yoh 


where R, is the fixed position of the tank center. Since y, is 
constant, to first order ; 


R; ~il (ir — Â cos Yo— jite 008 Yo sin yg + 
kl,(3,— 8) sin yo 


Again keeping only first-order terms leads to the following 
equation: 


oX R/ «ox R4 ili cos yo—Ji (bcos ye Y sin y) +klð sin yo 
Now combining the last two equations gives 
ox R/ +R sex R- ingly cos yo— Jut; kl; sin "Yo (5) 


This equation shows, as could be predicted physically since 
no viscous force is assumed between the tank and the fuel, 
that the airplane rotation affects only the motion of the tank 
center.. 

From equations (4) and (5) the necessary translational 
velocities can be obtained for the translational kinetic ener- 
gies. The rotational velocity of the airplane is simply c. 
The spinning motion of the fuel about the vertical axis is 
ignored; then, the rotational velocity of the fuel may be 
given by the components 7 and č. The two corresponding 
horizontal axes of rotation through the fuel center of gravity 
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are principal axes of the spherical segment of fuel; conse- 
quently, no product-of-inertia terms occur in the fuel rota- 
tional energy. Also, since the airplane center of gravity is 
in the airplane symmetry plane, Jyz=Jyy=0 and only the 
Iz product of inertia will appear in the airplane rotational 
energy. 

By use of equation (4), the airplane velocity can be shown 
to be 


V, (U +u+6h,,—oR,) (o $R, — oR) + 


k(w4-óH, —óR,) (6a) 


When equation (5) is substituted into equation (4), the 
velocity of any particular fuel center of gravity is 


V,—iU +u E 6R,, VR, V ad; cos vo) - J(c E, —óR, — 


fds) +k(w+oR,,—6R,,+ al; sin yo) (6b) 


If V is the magnitude of the translational velocity of the 
center of gravity of a rigid body, J’ is its moment-of-inertia 
tensor about the axes through its center of gravity, and w 
is the rotational velocity of the rigid body, the kinetic 
energy is 

1 nail / z £ 
Er=5m V+ 5x wt Ly ey Ez w 


22 yz’ wyo,— 
2F yz’ oa, — 2I ey’ wro ` (7) 


Thus, for the kinetie energy of the airplane, substitution of 
equation (6a) into equation (7) gives 


E, = l(U+ut oR, VR, + (o cR. — 6R, )!4- 


G 


r! Ly: 


Lege 





(w+ oR, —6R, yp 


and, for the Kinetic energy of each fuel, substitution of 
equation (6b) into equation (7) gives 


E= (U +u +R, — VE, 4- ñl, cos Yo + 
(týR — 6R, id) (w+ oR, — 6H. + sin yo)’ + 
— RE (8b) 


k 
For the total kinetic energy, E,— E; 4-2 ; Ex; therefore, 
fio! 


equations (2), (3), and (8) may be used in equation (1) to 
obtain the (2k+6) equations of motion. However, it must 
again be recalled that the coordinate system is rotating. 
The whole system is therefore subject to an additional 
gyroscopic acceleration since the time-derivative operator 
contains an additional gyroscopie term (see reference 6) 








when the components of the velocity (or any vector) are 
taken in the rotating system: 


Bx 
uF-—r-LexXr 
õi T 


Thus the gyroscopic acceleration acting on the whole rotating 
system is 


eXr =0X (iU +0) = (eX DU =jyU— kiU 


The effect of this acceleration can be brought into the equa- 
tions of motion by considering the inertial reaction of the 
total mass to this acceleration i an additional applied force. 


If the total mass is m,=m,+ Sim, this reaction has the 
folowing components: 


Y'—-—m,Uy| 
(9) 


Z- mU ê 


In addition there is the inertial reaction torque M, on the 
fuel; this torque acts about the tank center and is caused by 
the acceleration of the tank center. For each tank, the 
vector reaction torque is 


Myf =R} — R)xm,U(—- Hy kô) 
=m;U [(—il, sin yo+ kl, cos y) X (—jý + kô) + 
Second-order terms] 
e mjl,U (ly cos yat- jê sin yot ky sin yo) 
Since M;'—(M/), cos yot (Mj), sin yo and M,'—(M,),, 
(102) 
(10b) 


Mj =m d, U} (cos*y,--sin*y) — mU 
M, —ml,U sin y, 


The forces and moments in equations (9) and (10) must 
be added to the weight and aerodynamic forces to obtain 
the Q: in equation (1). 

For convenience, the results of equations (2), (3), and (8) 
are as follows: 


1 £ 2 
E=3 921md G7 as) (118) 
E,—7* (Ur u-- iR oR, (o $e, — OR. Y + 
. . k ` 
(wt oR, — 8R: y] TÉ s KU wt 8R,— 
oR, + 1,1; COS yo + (o - VR. ,- OR, E 
‘ : CR alode 
(wt GR, — UR, id, sin vI} E dea 
Ig x Iz: FE 1. te od is 
S € Is): Uh) 
2 2 2fei ^! 
(11b) 
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Equations (9), (10), and (11) may now be used in equa- 
tion (1) to obtain the equations of motion. It should be 
noted diia 9... 9. DM and Dr 

ot O(U--w) òu Oy Ov O2 Ow 
ample, to obtain the equation of motion in the z-direction, 
note that 


For ex- 


OE, 8E, _ 


ru 








Then, equation (1) may be written as follows: 


dnas om 
dt\ ot / dt\ ou 
=m. (i+ OR. — 9I, ) - 2: mF OR. VR y+ 


3l; COS Yo) 
oy k 
=myi+5( ma, 2m, ) (m. 
fei 
k ko 
zin, jmd, cos Yor (12) 





The position of the origin of airplane coordinates has not 
yet been specified. Equation (12) and the similar equations 
obtained in the other degrees of freedom suggest that the 
position of the origin be determined by the following three 
conditions: 


k 
mR, d 25 m,R,,=0 (13a) 
{=l 
E 
malty, + 2 m,H,,—0 (13b) 
mR tZ mR, =0 (136) 
f=ì 


Equations (13) imply that the origin is at the position of the 
total center of gravity when the fuel mass is treated as being 
concentrated at the tank center. This choice of the origin 
greatly simplifies expressions such as equation (12). The 
physical reason for this choice is again that the fuel docs not 
rotate with the airplane; thus, a force acting on a line through 
this point, the center of gravity where the fuel reaction is 
assumed concentrated at the tank center, will produce no 
rotation of the airplane. 

The following substitutions will also greatly simplify the 
writing of the final equations of motion: 


Lm Toms [(R) RIEA m 0) 08 


Io Li deme (Ra) H (Ra) HX ml (Re) + (Re) 
(14a) 


k 
Il; + mel (Re) t (S) 22 m US) 0] 
Lom ID! m REA mR: Ra 
-I 


Note that the quantities defined by equations (14a) are the 
total moments and product of inertia about the origin of the 
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airplane coordinates when the fuel mass is assumed to be 
concentrated at the tank center. Finally, the necessary 
moments of inertia of each fuel pendulum about the tank 
center are 


Isl tmt | (14b) 


Famam 
Without loss of generality then, equations (13) and (14) 
are used in the equations of motion obtained by substituting 


equations (9), (10), and (11) into equation (1). The general 
equations of motion eau now be given as follows: 


k 

mata mjljij, cos yo X 

l k 

milth— U8) + 2s madir Sin y= Z 

=1 
DN k 7 
Iy 8 (nt, Re 23 m Tl, ó— 
fei i 


k 
(m.s tu tA ma) i 


<h (Jt, cos yo— Re, sin ym M p 059 
{z1 
Ld mh [gn Ht cos yot (h— U8) sin yo+ 
H,ó sin yot OR, cos Yo— Fz, sin yo) — 
R ý cos Yo] =0 
ere ae: Ave 
mí EU) 23 mj, 
^ - k . 
Tey Lao (ns, Regt 2) mR, Baji = 
: * at 
zi ml (Rt e uh, cos yo) =N 
Id Ici - (nes 23 m) | sy 


k 
21 mij, + Roi, sin Yo) =L 


I, hd ml(gti—i—UP-- Rab Rao) =0 


In equations (15a) and (15b) only the fuel equations for the 
first fuel tank have been written. In each set there are E 
similar fuel equations. As has been previously stated, the 
forces on the right-hand sides of these equations are the 
applied forces and the weight and aerodynamic forces. 
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SIMPLIFYING ASSUMPTIONS 


The equations of motion have been separated into what 
would generally be considered the longitudinal motions, 
equations (158), and the lateral motions, equations (15b). 
In the ordinary six-degree-of-freedom case, as can be shown 
from considerations of symmetry, no cross-coupling terms 
exist between these motions in the aerodynamic forces (see 
reference 7). Although such terms are known to exist in 
practice, they are small and generally neglected. However, 
many eross-coupling terms occur between equations (15a) 
and equations (15b) because of the fuel motions, even when 
the aerodynamic coupling forees are ignored. The magni- 
tudes of these fuel forces can be seen to depend on the masses 
of the fuels, the vector positions of the fuel tanks, the “pendu- 
lum length” (i. e., the radius of the tank and the height of 
the fuel), and the accelerations involved in the fuel motion. 
Increasing the magnitude of any of these parameters will 
increase the effect of fuel motion. For most airplanes these 
fuel forces will be relatively small, but the present investiga- 
tion is primarily concerned with all first-order fuel effects. 
In any particular case, the actual magnitudes of these forces 
may of course be obtained by inserting the values of the 
previously mentioned parameters. 

Symmetrical fuel distribution.—Some of the terms in 
equations (15a)—for example, the œ and $ terms in the 8 
equation—are essentially product-of-inertia terms arising 
from unsymmetrical distribution of the fuel about the y—0 
plane (i. e., the plane of symmetry). That this is so can be 
seen if the fuel tanks are assumed to be distributed symmetri- 
cally with respect to the symmetry plane; that is 


k k k 
2; mR, Rh, =>) mR R,—» m,R, —0 (16a) 
j=1 fo si ON Fat d 

and from equations (13) 


R, —0 (16b) 


Therefore, these terms vanish for symmetrical fuel distribu- 
tions. In most cases the fuel will be symmetrically distrib- 
uted, and substitution of equations (16) into equations (15) 
yields the equations of motion for symmetrical fuel distri- 
hution: 


k ` 
mer Mel pity cos Yo X 
f= 


k 
mlw—Cb+t >) Mie SM yo—Z 
f-i 


ay ko . i 
Ir d 0 UL, cos yo— Fz, sin yom li, — M 
fz1 ; f 
| . ; (7a) 
[hcc mil[goid cos yot (é— UB) sin yo+ 


Fi, ó sin yot 6 (1, cos yo— E. sin yg)— 
Py cos yo — 0 


Tí2182—34 31 











3 k * 
mor Ep) mdt =Y 
=i 
£ 2k S 
Inb—Ixyzó E MAR, trt R, ñr cos y) =N 


E ad k B 
Ixó— Ted + 23 mA. Cpt Ry iy siny)=L (7b) 


I fit mh(gti—£—Utbt--Rj— R.é) =0 


In equations (17) even though the terms arising from unsym- 
metrical fuel distribution have vanished, some cross-coupling 
terms still remain between equations (17a) and (17b). These 
terms occur in the y equations of set (17a) and in the ¢ and V 
equations of set (17b). The significance of these terms is 
evident since each contains a factor R, Thus, these terms 
arise when the airplane has fuel tanks with centers not in 
the plane of symmetry, even though they are symmetrically 
distributed with respect to this plane. For example, they 
would arise for wing-tip tanks. Physically, these terms 
clearly give the interaction between the longitudinal fuel 
motion s, in the wing-tip tanks and the airplane rotation 
&bout the vertical axis, which consists of the lateral motions 
$ and V. For example, assume for simplicity that ya=0; 
then, à yawing acceleration of the airplane will cause a lon- 
gitudinal fuel acceleration 7, in the wing-tip tanks, and vice 
versa. 

From this discussion the 7 motion appears to couple the 
lateral and longitudinal airplane motions even for the per- 
fectly symmetrical fuel distribution described by equation 
(16a). However, the fact that this coupling does not occur 
can be seen by considering any pair of symmetrically placed 
and loaded tanks. Designate the 7 motion in this pair of 
tanks by m and qw. Then, the y equations in equations (17a) 
show that a longitudinal horizontal acceleration of the tank 
gives rise as expected to ņ accelerations. Since the system 
is linear, this portion of the n motion may be considered 
independently, and because of the symmetry of the two 
tanks it is seen that 7,=m for the portion of the y motion 
arising from the longitudinal motion. Therefore, in the V 
and é equations of set (17b) the effects of this 7 will vanish 
since E, — —H,. In a similar manner the laterally caused 
7 motion can be shown to have no effect on the longitudinal 
motion. Essentially the argument is that the 4 motion for 
each pair of tanks can be split up for perfectly symmetrical 
fuel distribution into symmetrical and antisymmetrical 
motions. The symmetrical portion of the 7 motion for each 
pair of tanks couples only with the longitudinal motion; the 
antisymmetrical 7 motion couples only with the lateral 
motion. Thus, for perfectly symmetrical fuel distributions 
the n equations of set (17a) could be combined with set (17b), 
only lateral degrees of freedom in the 7 equations being used 
(since the symmetrical portion of the y motion is of no 
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interest); or the 7 equations may be used as shown with 
equations (17a), the lateral degrees of freedom in the y 
equations being ignored. 

Fuel tanks centered in symmetry plane.—Many airplanes 
have large internal fuel tanks which are centered in the 
airplane symmetry plane. For such airplanes the equations 
of motion may always be separated into independent lateral 
and longitudinal modes since R,,=0. Using this value in 
equations (17) gives the equations of motion for tanks 
centered in the symmetry plane: 


k 
müt S mily cos yo=X 
i 
AS $2 k * 
m,(w—U8)+ >) miii, sin You Z 
f=1 


7 k . 
Ir$t2i (Ra, cos Yo— Fs, sun yg mjlgm,-— M 


(18a) 
I, mil [gm + & cos yo+ (b — U6) sin yot 
(Ra cos Yo— Fi, sin y3:)8]—0 
: "EN: PP 
m(v TUD—2 mj, f= 
Ip—1x;9—25 mj,R,p—N 
(18b) 
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Each set of equations, (18a) and (18b), contains (k4-3) 
variables; the two sets can be seen to be independent of each 
other since the ¢ motion of each tank couples with only the 
lateral motion and the 7 motion couples with only the longi- 
tudinal motion. 

For the case of a single fuel tank at the airplane center of 
gravity the modification of equations (18) is obvious. Then 
R,=R,=0, and all coupling between the rotational motion 
and the fuel motion vanishes; that is, all the fuel terms in the 
rotational equations vanish and all rotational terms in the 
fuel equations vanish. For an aspherical tank the rotational 
coupling in this case will be small. 


LIMITATIONS INHERENT IN THE APPROXIMATIONS 


Before proceeding to the application of equations (18b) 
it is appropriate to consider somewhat more explicitly the 





assumptions involved in the indiscriminate dropping of all 
second-order terms which appeared during the derivation 
of the equations of motion. In this connection the correc- 
tion, arising from the airplane accelerations, to the constant. 
acceleration field g involved in the pendulum potential en- 
ergy should be considered. The assump tion which is implied 
in neglecting these accelerations is that the accelerations of 
the tank centers are small with respect to g. 

If this and previous approximations are considered, it 
can be seen that three essential assumptions were made in 
dropping second-order terms: 

(1) The fuel and airplane angular displacement variables 
are small enough so that the angle approximates its sine. 
However, this approximation sometimes took the form that 
the angle was much Jess than 1 radian. 

(2) The disturbance velocities are much less than U, and 
products of the linear or angular velocities can be ignored. 

(3) The accelerations of the tank centers must be small 
compared with g. i 

Strictly speaking then, the statement that the equations 
of motion (15), and also the simplified equations, are ac- 
curate equations of motion to first order is to be taken to 
mean that the motions to which these equations apply are 
restrieted by the preceding three conditions. Thus, the 
equations would appear to remain accurate at least at the 
beginning of a disturbance. Moreover, when the motion 
becomes large enough that these assumptions break down, 
the fundamental physical assumption that the fuel may be 
considered to move as a rigid body also breaks down; there- 
fore, nothing can essentially be gained by keeping higher- 
order terms in the mathematical expressions. 

Since the pendulum motion is little changed even up to 
angles of 30° to 40°, it could be expected that aside from 
splashing effects these equations should remain a good to 
fair approximation even at such angles. On the other hand, 
even the splashing effects, although they would introduce 
some damping and change the inertial characteristics of the 
pendulum somewhat, could certainly not be expeeted to 
cause the general assumptions to break down completely 
for fuel motions up to angles of 30° to 40°. Therefore, the 
equations of motion derived are assumed to present a fair 
picture of the disturbance motion even up to fuel displace- 
menís of tbis magnitude. 


APPLICATIONS TO SEVERAL CASES 


NONDIMENSIONAL EQUATIONS FOR TANKS IN SYMMETRY PLANE 


The equations of motion (18b) have been applied to the 
lateral motion in several cases with two fuel tanks in the 
plane of symmetry. In these cases the lateral motion can 
be considered independently. The applied forces are the 
weight, the usual aerodynamic forces linear in the disturb- 
anee velocities, and any disturbing forces that may be 
present. In order to put the equations in nondimensional 
form, the nondimensional lateral airplane equations are used 
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as obtained in reference 8. "The fuel equations are made 
nondimensional by making the standard transformation to 
uondimensional time derivative, as in the airplane equations, 


m4,U* 
b 


and then by dividing through by —“=—. The resulting 


nondimensional expressions in the following equations are 
defined in equations (20): 
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The derivatives er. and Cr, were assumed to be zero. 


METHODS OF SOLUTION 


In the present case the two fuel degrees of freedom intro- 
duce two additional oscillatory modes into the characteristic 
solution, in addition to modifying the original airplane mode. 
The motion will therefore be a combination of three oscilla- 
tions (aside from the less important exponential modes), but 
just knowing the three oscillatory roots is insufficient to indi- 
cate the type of motion since the relative magnitudes of the 
oscillatory modes must also be known. For this reason 
motions must be caleulated in order to see the actual effects 
of the fuel motion. However, in several cases the charac- 
teristic roots were also found in order to facilitate the inter- 
pretation of the motions. These cases will be discussed 
subsequently. 





The most convenient method for the analytical solution 
of a set of linear ordinary differential equations such as 
equations (19) is probably the Laplace transform method 
(see reference 8). However this method is extremely eum- 
bersome and difficult to check since it involves the expansion 


` of fifth-order determinants in which the elements are often 


quadratie functions of the characteristic root. Therefore, it 
seemed preferable to use some step-by-step method which 
would be more amenable to machine computation. 

Reference 9 gives a matrix method for getting the step-by- 
step solution of a set of linear ordinary differential equations. 
When applied in the present case to equations (19) this 
method results in a simultaneous solution for the motion in 
each of the five degrees of freedom and also for the motion 
in Dé, Dy, Dj, and Dg; The calculations were carried out 
on the Bell Telephone Laboratories X-66744 relay computer 
in use at the Langley Laboratory. The essential details of 
the method are given in appendix A, 


SOLUTIONS FOR SEVERAL CASES 


The two basic cases for which motions were calculated 
Were case A, a present-day high-speed airplane with two fuel 
tanks satisfying the conditions for equations (19), and case 
B, which corresponds essentially to case B of the model used 
in reference 2. The essential parameters for these two cases 
are given in tables I and II. Table II gives the conditions 
for case A when both tanks are one-half full (Ar) and when 
the fuel height equals one-half the radius (À3) and for case 
B when the fuel heights in both tanks are 2 inches (B5), 3 
inches (Bi), and 4 inches (B,). 

In case A the tanks are spherical, somewhat over 4 feet in 
diameter, and centered on the body axis approximately 4 
feet in front of and behind the airplane center of gravity. 
The flight conditions are given in table I. The fuel weight 
in the half-full condition is &pproximately 25 percent of the 
total weight. 

In ease B the tanks are spherical, centered in the plene of 
symmetry slightly less than 5 inches below the airplane axis 


TABLE I 
STABILITY DERIVATIVES AND MASS AND GEOMETRICAL 
CHARACTERISTICS FOR TWO CASES CONSIDERED 











i Parameters | Case A Case B f 
j : 
Cr,.. AMT DD —0.80 
Mp TREE | oi | *0.I7 
PO a aaas aan EN telo eti el V DA LL —0. 126 —0.14 
E DO. es ----| 0.01552 5—0.05 | 
Ct —0. 342 —0.30 
—0.28 —0. 16 
0. 0796 0.30 
6, 910 11.25 
130 2.07 
| 28 40 
0 -—1i 
0. 00138 0. 002378 | 
| 704 6L 5 to 74.25 i 
: 4d 














* In case Bis, Cag=0.29. 
b Actually, slightly different values of Ca, were used for each of the subcases of case B. 
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and 4 inches in front of and behind the model center of 
gravity. The diameters are 8 inches and the total fuel 
weight in the half-full condition is approximately 46 percent 
of the total weight. The flight conditions are the same as 
in reference 2. 

Motions were calculated for certain subcases of the basic 
cases which were obtained by varying certain significant 
parameters. By comparing the resulting motions an attempt 
was made to evaluate the effect of varying such factors as 
amount of fuel, position of tanks, and relative natural fuel 
and airplane frequencies on the disturbance of the airplane 
motion caused by fuel motion. Also various initial condi- 
tions were considered to show the effect of initial conditions 


TABLE II 
FUEL DEPENDENT PARAMETERS 






































Parameters Ay Aa Bs Bı Bi 

Radius of forward 

tank, It ees coheed a= 2.15 2.15 0. 333 0. 333 0.333 
Radius of rear tank, ft. 2.12 2.12 0. 333 0. 333 0. 333 
Fuel height in for- 

ward tank, ft-------- 2.15 1.075 0. 1667 0. 250 0.333 | 
Fuel height in rear 

fink ftorni 2.12 1.00 0. 1657 0. 250 0. 333 
Rip ft...----- 4 » 3.5 0. 333 0.333 0. 333 
Ray, {t..--~ —4.1 —0. 333 —0. 333 —0. 333 
Ray ft....- 0.0123 0. 407 0. 407 0.407 
Ruy, (t.------------ nee —0. 0144 0. 407 0. 407 0. 407 
Weight of forward 

uel, Ib....---------- I, 480 462 151 3.07 4.85 
Weight of rear fuel, Ib. 1, 068 334 1.51 3.07 4. 85 
Total weight of air- 

plane and fuel, 1b. - . 9, 518 7, 766 14. 27 17. 39 20. 95 

E S EEEE 0. 19 0. 19 1. 062 1.178 1.174 
Ix, slug-ft?..- x 1, 360 1,300 0. 1081 0. 1243 0.1428 
Iz, slug-ft?. ....-------- 7, 340 7, 708 0. 2098 0. 2206 0. 233 
Izz, slug-ft?. z 262 | 277 0 0 0 
ly ft... K 0. 806 1.45 0. 225 0.174 0.125 
Ce eee » 0. 795 1.43 Q. 225 0. 174 0.125 
hi, slug-ft?_ zi 85.2 40.0 0. 00315 0. 00507 0. 00668 
Ayslugtt soe 59.7 28.1 0. 00315 0. 00507 0. 00668 
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motion scems to have more damping than the natural airplane 


on the resulting motion. In some cases an initial disturbance 
in sideslip was assumed, and in other cases an initial fuel 
disturbance was assumed. An initial sideslip of 5° and a 
fuel displacement of 10° were arbitrarily chosen as standard. 
Since the equations are linear, multiplying the initial dis- 
placements by a common factor simply multiplies the result- 
ing motions by the same factor. In this connection it must 
be emphasized that, if at any time the calculated motion in 
any degree of freedom becomes too large to satisfy the 
limitations previously diseussed, the following motion is 
meaningless. For example, if a 5? displacement in B gives 
rise to a fuel motion much greater than 30? to 40°, the accom- 
panying 8 motion is meaningless because the assumption of 
small displacements is violated. However, if multiplying the 
fuel motion by some arbitrary factor, for example 2/5, will 
bring its peaks down to less than 30° to 40°, then the 8 motion 
resulting from an initial £ disturbance of 2° can be obtained 
by simply multiplying the previous 8 motion by 2/5 also. 
The effect of large fuel displacements must be discussed 
qualitatively. 

The motion in sideslip and the motion of the two fuel 
pendulums in the various subcases are shown in figures 2 to 
12. Comments on these motions are presented to facilitate 
interpretation of the figures. The period of each fuel pendu- 
lum is called the natural fuel period. The period and damp- 
ing of the airplane, the fuel being disregarded, are called the 
natural airplane period and damping. 

Case A,— The natural fuel periods for case A; (half-full 
tank) are approximately 1.66 seconds and the natural period 
of the airplane alone is 1.40 seconds. Damping to half- 
amplitude occurs in 2 cycles. 

"The motion in figure 2 is for initial 8,—2^. The early 8 
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Fiqure 2. -~Present-day high-speed airplane with tanks one-half full (case Az). Sideslip and fuel motions following initial sideslip, 802 Z^. 
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mode. The disturbance arising from the fuel modes is evi- 
dent after 2 eycles. That the irregular residual oscillation 
of amplitude 4? to 5? is essentially due to the fuel modes 
is seen from the fact that the dominant period in the later 
motion is approximately 1.6 seconds. Notice that in this 
case a 5° initial fa would almost immediately cause fuel dis- 
placements of over 80?, so that the following motion would 
he radically changed. 

The motions shown in figure 3 for ({)o=(t2)o=10° are 
quite regular and indicate one dominant mode in each motion. 
The fuel period is 1.6 seconds. The airplane period starts 
at 1.4 seconds. builds up to 1.7 seconds, and averages 1.6 
seconds. The amplitude of the sideslip motion is very small. 
The largest such motion which could occur for this type of 
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disturbance would be for initial (¢:)o= (12)9—30? and would 
give 8 amplitude slightly more than 0.1°. 

The small amplitude of the sideslip motion in figure 3 was 
surprising. It was conjectured that for this fuel configura- 
tion the fuel displacement ({:)»=10° and ([3e— —109?, 
corresponding essentially to an initial yawing moment, might 
be more effective in inducing an airplane oscillation. (See 
fig. +.) Apparently, this configuration is more effective 
inasmuch as the sideslip motion now builds up to an ampli- 
tude of approximately 0.4°. The energy necessary to induce 
this considerable *snaking" type of oscillation seems to be 
obtained initially from the rear-tank motion, which is in the 
proper phase relation with the sideslip motion to feed energy 
into it at the start of the motion. 
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Fir RE 3.—Present-day high-speed airplane with tanks one-half full (ease A;). Sidestip and fuel motions following initial fuel disturbance, (f1)o=({s)e=10°. 
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The fuel periods in the regular motion are slightly over 


1.6 seconds. The airplane period increases from 1.5 to 1.8 
seconds and has an average period of 1.6 seconds. 


period of approximately 1.7 seconds. Both fuels start with 
a period of approximately 1.5 seconds, which increases to 
1.7 seconds. f l 

Case A;.— The natural fuel periods in case A, are approxi- In figure 6, (f)e——10? and (f)>=10°. The sideslip 

mately 1.52 seconds and the natural airplane period is 1.49 | builds up to a fairly regular oscillation of 0.4° amplitude 

seconds. Damping to half-amplitude occurs in 1.3 cycles. with the period increasing from 1.5 seconds to 1.7 seconds. 
The motion in figure 5 is for initial 8,—0.5?. The early | The fuel motion has a period somewhat under 1.7 seconds 

sideslip motion seems to be of greater damping than the | with amplitude quickly building up to the limits where 

natural airplane mode. ‘The residual airplane motion arising | splashing must become important. 

from the fuel modes sets in very quickly and is a regular 

unstable motion of very large relative amplitude with a 





Case B;— The sideslip motion shown for case B, appears 
to be a normal damped oscillation for the first 4 seconds 
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FIGURE 4.—Present-day high-speed airplane with tanks one-half full (case Ai). Sideslip and fuel motions following initial fuel disturbance, (f:)o=10°, (¢2)o=—10°. 
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(fig. 7, 8,—5?), but then the peaks show a slight irregularity 
instead of damping smoothly. The period of the early side- 
slip motion appears to be somewhat over 0.9 second and the 
motion damps to half-amplitude in less than 2 cycles. "This 
motion is very close to the undisturbed airplane model 
tperiod of 0.92 sec and damping constant of 1% cycles). 
The natural fuel periods are 0.61 second. The fuel motion 
is very irregular and obviously contains considerable amounts 
of at least two characteristic modes. 

Cases B, and B,.— The motion in cases B; and B, shown in 
figures 8 and 9, respectively, for 8,—5? is very much the 
same as in the previous one, except that the disturbance of 
the airplane mode in the sideslip motion appears somewhat 
more pronounced as the amount of fuel increases. 

Figure 10 (case By, (n)g— ({2)o>=10°) shows that the motion 
in sideslip resulting from the fuel displacement is much more 
irregular than in case A. The dominant mode corresponds 
to a fuel frequency, but apparently the airplane mode is 
present with considerable amplitude. The maximum oscilla- 
tions are approximately +%°. The sideslip motion in this 
case was much more irregular than for the corresponding 
initial conditions in case A. It was conjectured that this 
might be caused by the fact that in this model both tanks 
are below the X-axis, so that the coupling of the fuel motion 
with yawing and rolling motions does not have the same 
phase relationship as in case A where one tank is above and 
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one is below the X-axis. Therefore in case B4, shown in 
figure 11, the front fuel tank was assumed to be above the 
X-axis, all other conditions remaining as incase B,. In this 
case the general type of motion does seem to resemble that 
in figure 3. The sideslip, which builds up to +°, shows 
a snaking at the fuel frequency. 

In case A the fuel natural periods are very close to the 
airplane period. In case B, however, the fuel period is ap- 
proximately two-thirds of the airplane period. In case Ba, 
shown in figure 12 for B85—5*, the value of Cag of the model 
has been arbitrarily changed to give the model a period very 
close to the fuel period of 0.66 second. Comparison with 
figures 2 and 5 shows that the motion in this case is very 
much like the motion in case A. 

Transverse accelerations.—In evaluating pilots’ reactions 
to snaking oscillations, the magnitude of the transverse 
accelerations involved in the oscillations has been found to 
be an important factor. Acceleration amplitudes above 
0.025g are found to be bothersome, and amplitudes above 
0.08g are considered very unsatisfactory. Calculations of 
the transverse accelerations involved in several of the previ- 
ous motions were carried out. The magnitudes of the accel- 


eration peaks in the residual oscillations were found to be 
approximately 0.04g to 0.05g. The actual motions are not 
shown since all the airplane oscillations are essentially of the 
same type as the 8 motions. 
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FIGURE 5.—Present-day high-speed airplane with fuel heights one-half the tank radius (case As}. Sideslip and fuel mations following initial sideslip, £$ =0.5°. 
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} FIGURE 6.—Present-day high-speed airplane with fuel heights one-half the tank radius (ease Az). 


Sidcslip and fuel motions following initial fuel disturbance, (£))o —10^, (£3)ox 109. 
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Frerre 7.—Free-fiying airplane model with fuel heights one-half the tank radius (case Bz). Sideslip and fuel motions following initial sideslip, 80—5?. 


DISCUSSION OF RESULTS are given as follows in terms of periods and times to halve 
or double amplitudes, in seconds: 


The motions described in the preceding section are suffi- 
Case As, natural modes: 


cient to give a fair picture of the types of possible fuel effects. 
Moreover, since each motion is just a superposition of the 
characteristic modes of the total system, these motions are 
often easier to understand if the characteristic roots are P,=P,=1.52 
known. Physically it is clear that the characteristic modes 
will not differ much from the natural (uncoupled} modes Case Ay, total system: 
when the interaction between the airplane and fuel is small. ` 
Comparison of figures 2, 5, and 12 with figures 7, 8, and 9 P,—147 T4—61 
indicates that the interaction between airplane and fuel is 
strongest when the frequency of the airplane is elose to that 
of the fuel, as might be expected from comparison with the P " = 

Sq : s 2771.67 T,=4.56 
resonance phenomena exhibited by an oscillator driving a 
system at its natural frequency. For this reason the charac- 
teristic modes of the total system were calculated in cases 


P,—1.49 7,,—1.91 





Case Bia, natural modes: 


A. and Ba, where the frequency ratio between airplane and P,=0.66 T= 1.29 
fuel natural frequencies was practically unity. The natural 
modes are given for purpose of comparison. The results P,=P,=0.66 


272483 —54—32 
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FIGURE 8.—Free-flying airplane model with fuel heights three-fourths the tank radius (case B3). 


Case By, total system: 


D,—0.65 T,,—30 
P,—0.54 T,—1.5 
»=0.71 T=0.48 


In these cases of large interaction it is difficult to identify 
one of the characteristic modes as the airplane mode. The 
characteristic mode in which the period is changed least 
from the natural airplane period has been called the airplane 
characteristic mode. However, in figures 5 and 12 this 
mode is not obviously the dominant one, as the airplane 
mode is in figures 7, 8, and 9. 

The most important effects to be noted in these particular 
cases, where the fuel and airplane frequencies are equal, are 
that a characteristie mode which is very lightly damped with 
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Sideslip and fuel motions following initial sideslip, fex 5?. 


& frequeney close to the natural airplane frequency exists 
and that an unstable mode appears. In connection with the 
first of these effects, it would seem that, theoretically, 
certain initial conditions might be found that would excite 
mainly this lightly damped mode in the characteristie solu- 
tion for £, so that the resulting motion would be a typical 
snaking. Of course the required initial conditions might or 
might not be practical ones. 

The total characteristic modes were also calculated for 
cases B, and Baa to investigate the changes in motion caused 
by a hypothetical shift of one of the fuel tanks. The natural 
modes and the characteristic modes of the total system are 
given as follows for comparison: 

Cases B, and Bia, natural modes: 


P,— 0.84 
P,=P,=0.66 


Ty=1.64 
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Fievre 9.—Free-flying airplane model with tanks one-half ful] (case Ba). Sideslip and fuel motions following initial sideslip, 8e—3?. 


Case B, total system: 


P,=0.88 T4 —1.91 the part of the motion shown. The 8 motion in figure 11 

P,—0.63 TBAB seems to show the effeet of the unstable mode. It appears 

m in this ease that, when the fuel tanks are in front of and 

P,—0.53 T4—1.61 behind the center of gravity, the configuration with one tank 

(tane Bo vat ersten: above and one below the X-axis gives rise to an unstable 

Bi » i mode; whereas the configuration with both tanks below the 
P,=0.86 Ty=1.14 X-axis makes both fuel modes stable. 

P,—0.61 T.—19.5 A comparison of figures 2 to 4 or figures 5 and 6 clearly 

: 3 E shows that the initial conditions can have a very important 

P,=0.60 T,=10.2 effect, since the least stable mode does not necessarily become 


It is interesting to note that both fuels in figure 10 and 
also in figure 11 seem to follow the more stable fuel mode in 
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Figure 11.—Free-flying airplane model with tanks one-haif full and front tank eee, is be sured above the X-axis (case Bie). 
fno= (Polo 10°. 

had practically died out, the remaining motion might be 

considered to be the type caused by a fuel displacement. 


From this point of view figures 4, 6, 10, and 11 seem to 


dominant for a long time. The 8 disturbance was chosen as 
On the other hand the dis- 


be completely artificial. However, these motions are be- 
lieved to give a rough idea of the residual oscillations caused 
by fuel motion, at least insofar as magnitude is concerned, 
since, if the fuels were still displaced after the airplane motion 


a typical airplane disturbance. 
turbance of the fuels as an initial condition would seem to 


indicate that residual oscillations of the order of magnitude 
Actually, 


of M? fo M? might be expected in these cases. 
figures 10 and 11 would show oscillations of the order of 


€? for 30? fuel displacements. 
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FiGURE 12.—Free-flying airplane model with tanks one-half full and Cag adjusted to make natural fuel and airplane frequencies equal (case Ba). Sideslip and fuel motions following ~ 


initial sideslip, f= 5. 


Because of the lengthiness of the calculations, only eases 
Ay, Ag, and By, were carried out far enough to show the 
residual oscillations following a 8 disturbance. It is evident 
that the motion in case By, (fig. 12) resembles the motion in 
ease A (figs. 2 and 5) much more than it resembles the un- 
modified case B motion (figs. 7 to 9). The reason for the 
smaller relative fuel motion in case By, is probably the fact 
that the relative fuel mass is considerably larger than in 
ease A. In ease A, (fig. 5) the residual oscillation dominates 
the motion almost immediately. These results show that the 
importance of the residual oscillation depends mainly on the 
closeness of the natural airplane and fuel frequencies, that is, 
on the parameter which might be called the frequency ratio. 
Inasmuch as the previous discussion of the characteristic 
modes indicated that the frequency ratio was also the most 
important factor affecting the characteristic modes of the 
system, the frequency ratio generally can be seen to be the 
most important factor determining the disturbance of the 
normal airplane motion caused by the fuel. Moreover, case 
A indicates that for spherical tanks the fuel frequency may 


easily be of the same order of magnitude as the airplane 
frequency. Reference 4 indicates that the same is true for 
rectangular tanks and for arbitrarily shaped tanks of reason- 
able dimensions. Thus, even though the residual oscillations 
might occur at fuel frequencies, these frequencies would not 
be distinguishable from the norma! airplane frequency in the 
cases where the fuel effect is most pronounced, since in these 
cases the frequency ratio approaches unity. 

The effects of unstable modes cannot be understood with- 
out considering the nonlinear effects due to splashing of the 
fuel. For linear systems the presence of an unstable mode 
would imply thet the total system is unstable. It has been 
shown, however, that in an actual motion if the coefficient 
of the unstable mode in the solution for the airplane motion 
is very small compared with the coefficient of one of the 
stable modes, then the unstable mode will not appear in the 
early part of the motion. Now even in cases where the inter- 
action is weak, one of the fuel modes (with no natural damp- 
img assumed) may be unstable. In such cases the unstable 
mode in the airplane motions will be relatively very small, 
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while the unstable mode may be dominant in one of the fuel 
motions. Then this fuel motion may become very large 
before the effect on the airplane motion can be observed. 
At this point splashing will set in and the equations of motion 
no longer hold. Actually the energy lost in splashing may 
damp out the fuel motion so that the effect of this fuel mo- 
tion, and therefore of the unstable mode, may never ap- 
preciably appear in the airplane motion. 

In ease As, however, the interaction is large, and figure 5 
shows that theunstable mode doessoon become very important 
in the 8 motion. Here again the effects of large-amplitude 
fuel motion and, in particular, the nonlinear effects of fuel 
splashing should be considered qualitatively. In cases A; 
and As, because of the strong interaction, a sideslip of 5° 
would very quickly give rise to fuel displacements of the 
order of magnitude of 90° or greater. Physically, it is clear 
that the airplane motion is feeding energy into the fuel 
motion in this early part of the motion. Because of the 
large amplitudes of the fuel motion, a large component of 
this motion is in the vertical direction, so that a considerable 
part of the energy in the fuel motion will go into creating 
longitudinal airplane motion. Because of the symmetry con- 
ditions previously eited this energy will not be fed back into 
lateral motion, so that the fuel tends to stabilize the lateral 
motion by feeding some of the energy from the lateral into 
the longitudinal motion. Also the turbulence due to splash- 
ing will absorb energy which will then be lost altogether 
from the motion. 

One conclusion which can be drawn from this discussion 
is that the motion following large disturbances may be more 
stable than that following small disturbances. For example, 
the initial disturbances in figures 2 and 5 have been adjusted 
so that during the motion shown there is little splashing, and 
considerable residual oscillations are shown. But if the ini- 
tial disturbanees in these cases had been 5? or more in f, 
the energy lost in splashing in the early motion would pos- 
sibly be so great that the residual motion in 8 would be 
smaller than that shown here. This might explain why some 
airplanes which definitely showed troublesome fuel oscilla- 
tions were reported to be more stable in conditions of large 
atmospheric turbulence than in slightly turbulent atmos- 
pherie conditions. 

A more important conclusion is that the fuel can be used 
to stabilize the airplane motion by introducing turbulence— 
for instance, by use of appropriate baffles. This can be seen 
by noticing that the early part of the 8 motion in figures 2 
and 5 is very stable. But if most of the energy fed into the 
fuel in this part of the motion were converted into heat 
through turbulence, then as has been pointed out this energy 
could not be fed back into the airplane motion and the 





REPORT 1080 -NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


residual oscillations would not appear. Since the amount 
of energy lost in turbulence cannot be calculated analytically, 
it would seem that an experimental investigation of the 
effects of honeycomb or other turbulence-inducing baffles 
on the airplane stability would be desirable, especially in 
cases where the airplane and fuel natural frequencies are 
approximately equal. 

Finally it is possible in a strongly unstable case of residual 
oscillations, such as shown in figure 5, that the fuel may 
lose just enough of its energy in splashing to reduce its 
amplitude to where the motion is again smooth. Then, 
because of the instability of the system for small motions, 
the amplitude might again begin to build up. In this way 
continued oscillations of a more or less regular nature would 
occur when the calculations neglecting splashing show un- 
stable motion. This result is important because it shows 
that somewhat irregular small-amplitude oscillations can be 
expected when the ratio of the airplane natural frequency 
to the fuel natural frequency approaches unity, even for 
moderate fuel masses of the order of one-tenth the total 
mass or less. 


EFFECTS OF ASPHERICAL TANKS 


The calculations have been carried out for rigid-body 
motion in spherical tanks only. Actually this assumes that 
for small oscillations the fundamental wave motion in 
spherical tanks approximates rigid-body motion. This 
approximation only applies when the tanks are one-half full 
or less. This restriction is not too serious, however, since 
the fuel motion will generally have its greatest effect in this 
range. 

It is important to note that the potential energy of the 
fuels is simply the potential energy of a set of harmonic 
oscillators located at the positions R, Thus, the same 
general analysis will apply whenever the fuel motion in the 
tank can be represented in terms of harmonic oscillators 
with given effective mass and spring constant. Reference 4 
has already been mentioned as obtaining such a representa- 
tion for the fundamental mode of a rectangular tank. 
Usually the fundamental mode will be the most important 
and will involve the greater effective mass. It is conceivable 
that for long tanks the second mode might be of a frequeney 
closer to that of the airplane and in that case might be more 
important. In such a case each mode might be represented 
by a separate oscillator. As has been pointed out, the 
damping is mainly due to turbulence and will be more im- 
portant for aspherical tanks. For small motions, however, 
the damping may still be neglected. 

From the general derivation of the equations of motion, 
the most important result was the effect of fuel distribution 
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on the coupling of lateral and longitudinal motions. It is 
plain that these results apply strictly only to spherical tanks. 
Consider for example a tank of triangular plan form located 
in the symmetry plane and oriented symmetrically with 
respect to this plane. Because of the symmetrical orienta- 
tion it can be seen that, although sideslip motion will give 
rise to forward and rearward forces (because the pressure 
forces are normal to the diagonal surfaces), the forward and 
rearward motion will give antisymmetrical lateral forces 
which will cancel. In this case the result would be to feed 
energy from the lateral into the longitudinal motion; this 
condition would be favorable since the longitudinal motion 
is generally well-damped. For an unsymmetrically oriented 
tank of this type, energy could be fed back again from the 
longitudinal to the lateral motion and the problem would be 
quite complicated. 

In general, the results on coupling for spherical tanks 
would be valid to first order for such symmetrical plan forms 
as the rectangular or the diamond-shape ones. For any 
simple symmetrically oriented shape in the plane of sym- 
metry, & loss of energy from the lateral to the longitudinal 
motion might occur. This condition would be favorable. 
Finally, for tanks outside the plane of symmetry the same 
considerations would be valid if the tanks were symmetrically 
placed with respect to the symmetry plane and symmetrically 
shaped with respect to the plane through the tank center 
parallel to the symmetry plane. 








CONCLUSIONS 


The following conclusions may be drawn from the the- 
oretical analysis presented: 

1. Considerable disturbances of the normal airplane motion 
can be caused by fuel motion. 


(a) The most important factor determining the effect 
is the ratio between fuel and airplane frequencies. When 
these are equal, even moderate amounts of fuel (one- 
tenth the total mass or less) may cause considerable 
disturbances. 

(b) The most usual type of disturbed motion is a some- 
what irregular small-amplitude oscillation and the type 
of motion is strongly dependent on the initial conditions. 

(c) The effects of splashing will be to make the motion 
more stable, and the loss of energy in fuel turbulence may 
make it possible to increase the stability by artificially 
introducing turbulence in the fuel. 

2. The fuel motion may cause coupling between lateral 
and longitudinal motions. 

3. The derivation of the equations of motion for spherical 
tanks may be applied to any tanks where the fuel motion 
may be represented in terms of harmonic oscillators. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
LawGLEY Fiero, Va., November 21, 1950. 








APPENDIX A 


A STEP-BY-STEP SOLUTION OF THE EQUATIONS OF MOTION IN MATRIX NOTATION 


A matrix method for solving the equations of motion is 
given in reference 9. The first step in this method is to 
reduce equations (19) to a set of first-order equations by 
introducing as new variables Dó, Dy, Di; and Df, as de- 
scribed in reference 9 for Dé and Dy. This transforms the 
equations of motion into a set of nine linear first-order 
equations in the nine variables. In matrix notation the 
equations may then be written as follows when there are no 
applied forces Cy, Cn, or C: 


A(Dq) 3- Bg—0 (A1) 


where A and B are ninth-order square matrices and g is the 
column matrix (or vector), the elements (components) of 
which are the nine variables. In partitioned form, 


Lai Oas, 
p Fen s 


ege | 
Bz| n4 
Cii Eg 


q= latano, B Dra Dt Do, D] 


where I, is the identity matrix of fourth order; 0, and 05 
are zero matrices of order (4X5) and (5X4), respectively; 
and 
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Equation (A1) may be solved for Dg by multiplying 
through by the inverse of A: 
Dq— —A^BgzPq (A2) 
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Equation (A2) shows that, when there are no applied forces, 

the differentiation operator with respect to nondimensional 

time may be replaced by premultiplieation with the matrix 
z—AÀA-B. For example, note that 


D'a—D(Dg) - P(Dq) — P'q 


and thatsimilar relations would result for higher powers. 
Now, by Taylor's expansion, the value of q(s-- As) may 
be obtained from the value of g(s) by the series 


A 
q(sT-As)— TEE S Do) d D'q(s).. 
By use of equation (A2), this equation can be written 


(As)? p, 
vocas P Pt... |e) aa 
which is the fundamental recurrence relation used in the 
step-by-step calculation. In equation (A3), J, is the ninth- 
order identity matrix. The set of initial disturbances go 
being given, the magnitude of the step As will determine 
the number of powers necessary to obtain a given accuracy 
in the solution. Because of the relatively high fuel fre- 
quencies and because it was desired to obtain the motion to a 
rather large number of periods with reasonable accuracy, 
the series in equation (A2) was used to the sixth power with 
As approximately 1/20 of the airplane period. Thus, the 
matrix relation (À3) was 


g(s -As)— Qq(s) (A4) 
where 


garter st PAG py. +e P 
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